Classical Theory of High Field Atomic Ionization Using Elliptical Polarization 
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Important information about strong-field atomic or molecular ionization can be missed when using 
linearly polarized laser fields. The field strength at which an electron was ionized, or the time during 
a pulse of the ionization event are examples of such missing information. In treating single, double, 
and triple ionization events we show that information of this kind is made readily available by use 
of elliptical polarization. 

PACS numbers: 32.80.Rm, 32.60.+i 



I. INTRODUCTION 



The interaction between intense laser fields and gas 
phase atoms or molecules has attracted attention for the 
past two decades [i[ . Many strong-field phenomena have 
been observed and controlled, such as ionization of atoms 
[4I , dissociation of molecules [1] , generation of high har- 
monics [41 and creation of attosecond pulses |ff] . 

To make a theoretical description of the interaction be- 
tween an intense laser field and a multiple-electron atom 
is not an easy task. First, no analytical quantum mechan- 
ical solutions can be expected. Second, the laser electric 
field strength is comparable to the atomic Coulomb elec- 
tric field strength felt by a valence electron, so neither 
the laser field nor the Coulomb field can be regarded 
as a small perturbation and familiar perturbation the- 
ories cannot be used. Third, full-dimensional numeri- 
cal calculation of the time-dependent Schrodinger equa- 
tion (TDSE) is extremely demanding in computational 
resources and is effectively limited to the helium atom 

Simplified theories and models are thus desirable and 
have been developed to match rapidly emerging exper- 
imental results (for a review of such theories, see [M])- 
Among these theories, a semiclassical three-step recolli- 
sion model has been widely used to understand various 
strong-field phenomena heuristically [9|, |lO[. However, 
interesting questions (e^., the proper understanding of 
electron release times [ll[) are being raised in high-field 
atomic photoionization that have little or nothing to do 
with recoUision. Here we report theoretical calculations 
concerning single, double and triple ionization for laser 
intensities in the PW/cm^ range with elliptically polar- 
ized pulses and without recoUision. 

Two theoretical approaches have been used. An an- 
alytical theory is first developed for the sake of physi- 
cal clarity by extending the Simpleman theory 12] to 
include elliptical polarization. An illustration of Sim- 
pleman electron trajectories in an elliptically polarized 
laser field is shown in Fig. [TJ Numerical calculations 
are then performed using a classical ensemble method 
[l3[, which treats the entire system, both the laser field 
and the atom, purely classically but going beyond the 
Simpleman approach. We do this by including all forces 




FIG. 1: An illustration of Simpleman electron trajectories in 
an elliptically polarized laser field with ellipticity 0.5. Three 
close trajectories are shown, corresponding to three slightly 
different initial velocities Vy at the time of emission. 



(electron-ion and electron-laser) in obtaining solutions of 
the time dependent Newton equations (TDNE). 

The classical ensemble method has provided valuable 
insights into ionization dynamics |14| - [l9j . It takes into 
account naturally the occurrence or the non-occurrence 
of recoUision, so it can be used to study strong-field ques- 
tions with or without recoUision. Recently it has been ex- 
tended to include elliptical polarization and physical pro- 
cesses without recoUision [20] and good agreement with 
experiment [2l| has been achieved. 

The purpose of this paper is to show that elliptical 
polarization has the ability to uncover ionization infor- 
mation that is otherwise unreachable with linear polar- 
ization. For example, under what field strength was an 
electron ionized? At what time during the pulse was an 
electron ionized? Note that an electron may not be ion- 
ized at the pulse peak if the peak intensity of the pulse is 
higher than the intensity needed to ionize this electron. 
This is especially notable for intense laser fields that are 
capable of ionizing more than one electron. This ioniza- 
tion information cannot be directly and easily obtained 
with linear polarization. However, as we show, this kind 
of information can be straightforwardly obtained from 
the end-of-pulse ion momentum distributions obtained 
under elliptical polarization. Ion momentum distribu- 
tions can be measured using the COLTRIMS (cold tar- 
get recoil ion momentum spectroscopy) technique j22| . 



We will show that the momentum distribution of a singly 
charged ion reveals the ionization field and the ionization 
time of the emitted electron. The momentum distribu- 
tion of a doubly charged ion reveals the ionization fields 
and the ionization times of both emitted electrons. And 
the momentum distribution of a triply charged ion tells 
the ionization fields and the ionization times of all three 
emitted electrons. 

Such ionization information may reveal new ionization 
dynamics. For example, will electrons really follow the 
ionization fields and the ionization times predicted by 
independent-electron tunneling formulas [23|,|2J]? A re- 
cent experiment performed on argon has given a prelim- 
inary answer no, although the detailed physics is still 
under investigation |11| . 

The following questions can be answered by future ex- 
periments using elliptical polarization: How good are the 
commonly used tunneling formulas? What is exactly 
the role of the remaining electrons during an ionization 
process? Does this role depend on atomic or molecular 
species? Is this role the same for the first ionization pro- 
cess and for the second ionization process, and even for 
the third ionization process? 

This paper will be organized as follows. In section II 
we will first review the Simpleman analytical theory to 
explain how one can use elliptical polarization to obtain 
the above-mentioned ionization information. In section 

III we go beyond Simpleman theory and use the clas- 
sical ensemble method to perform numerical TDNE ex- 
periments and to test the accuracy of the analytical the- 
ory for an atom with three active electrons. In section 

IV results of the numerical experiments will be shown 
and compared to the predictions of the analytical theory. 
Summaries are presented in section V. 



II. SIMPLEMAN THEORY 

In this section we recall a simple analytical theory 
that links the experimental ion momentum distributions 
with the ionization information interested. This so-called 
"Simpleman" theory [12] has long been used to under- 
stand strong-field ionization processes, especially electron 
kinematics after emission from the parent ion. 

The Simpleman theory starts from the time that an 
electron is just emitted and neglects the ion core Coulomb 
potential. The ionized electron is treated as a classical 
particle and its motion is governed by classical mechanics, 
via solutions of the TDNE 



end of a pulse) is just 



F = md and F — qEL(t), 



(1) 



where EL(t) is the laser electric field. 

The momentum of the electron at the time of ionization 
is approximated to be zero: p{ti) = 0. The momentum 
of the electron at some later time t (for example, at the 



p{t) = / ait)dt 



EL{t)dt. 



(2) 



Atomic units are used. The charge of the electron is -1 
a.u. Note that for the intensities that we are interested 
in, relativistic effects are negligible and the magnetic part 
of the Lorentz force can be ignored. 

The Simpleman theory has mostly been used for single 
ionization with linear polarization. Here we will extend it 
to include elliptical polarization and to take into account 
single, double, and triple ionization. Let us start from the 
ionization of a single electron in an elliptically polarized 
laser field 

El (t) = Eof{t) [x sin(wi + (j)) + ye cos(wt + </))] , (3) 

with 



E^t) = Eof{t)sm{Ljt + 
Ey{t) = eEof{t)cos{ujt 



(4) 
(5) 



Where f{t) is the pulse envelope function, co the angular 
frequency, cj) the carrier envelope phase (CEP), and e the 
field ellipticity. We choose the x direction as the major 
polarization direction and the y direction as the minor 
polarization direction (recall Fig. [T]). 

Suppose an electron is ionized at time ti with zero ve- 
locity, and suppose that the ion core Coulomb attraction 
can be neglected after ionization. Then the momentum 
of this electron at the end of the pulse can be straightfor- 
wardly derived 



Pix = Ey{ti) w 0; 

ecu 

8 £ 

Ply = -E^{tl) = ±—El. 



U) 



LO 



(6) 

(7) 



An adiabatic condition has been applied. The dura- 
tion of the pulse is assumed to be much longer than one 
optical cycle. Then the end-of-pulse momentum of the 
electron does not depend on the detailed pulse shape. We 
see an interesting crossing relation: the end-of-pulse mo- 
mentum along the x direction {pix) is determined by the 
instantaneous laser field strength along the y direction at 
the time of ionization, Ey{ti), and the end-of-pulse mo- 
mentum along the y direction (piy) is determined by the 
instantaneous laser field strength along the x direction at 
the time of ionization, E,j; (ti ) . 

Note that for elliptical polarization, the peak field 
strength along the x direction is stronger than that along 
the y direction. Due to the fact that ionization proba- 
bility depends exponentially on field strength [23], ion- 
ization happens most probably around field maxima in 
the X direction. At such times, the field along the y di- 
rection is near zero. Therefore E^{ti) can be regarded 
as the instantaneous laser field at the time of ionization. 
Since Ex{ti) could be positive or negative, a "±" sign 
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FIG. 2: Schematic illustration of the ion momentum distribution along the y direction for single ionization (left), double 
ionization (center), and triple ionization (right). This figure is only for the purpose of demonstration, so it is not drawn to any 
scale and no further complications like the peak widths have been taken into account. 



has been added in front of Ei , which denotes the abso- 
lute laser field strength at the time of ionization and is 
always positive. 

We can see that the end-of-pulse momentum distri- 
bution of the electron, or equivalently of its parent ion, 
should have a single peak structure centered at zero along 
the X direction and a double peak structure centered at 
±eEi/uj along the y direction. Thus we have obtained a 
clear relation between the peak positions of ion momen- 
tum distribution along the y direction and the ionization 
field of the electron. 

If the ion momentum distribution is measured experi- 
mentally using COLTRIMS, projection of ion momentum 
onto the y direction should give two symmetric peaks, as 
illustrated in the left panel of Fig. [5] From the positions 
of the two peaks, noted as ±P with P > 0, one gets the 
ionization field of the electron as 



E, = ^P. 

e 



(8) 



This formula demonstrates why elliptical polarization has 
the ability to uncover ionization information unreachable 
with linear polarization: elliptical polarization provides 
an additional dimension, which contains information. 

The corresponding ionization time can be inferred as- 
suming a smooth pulse shape, e.g., a gaussian or a sine- 
squared pulse envelope. Figure [3] illustrates the idea us- 
ing a sine-squared pulse envelope and ellipticity value 0.5, 
which is the value used for illustration throughout this 
paper. Note that in general, there is no one-to-one corre- 
spondence between Ei and ti. The method that we have 




FIG. 3: Illustration determining ti after having obtained Ei. 
A sine square pulse with ellipticity value 0.5 is shown. The 
position of the first field peak corresponding to Ei is used as 
an estimate of ti. 



adopted to get ti in this paper is to use the position of 
the first field peak corresponding to Ei , as illustrated in 
Fig. [21 CEP phase will be averaged out. 

Next, let us move one step further for double ioniza- 
tion. Double ionization can be roughly divided into two 
categories, namely, sequential double ionization (SDI) 
and non-sequential double ionization (NSDI). SDI im- 
plies that the two electrons are ionized one by one with- 
out noticeable mutual correlations. NSDI means that the 
two electrons are ionized almost at the same time with 
substantial mutual correlations. Recollision is generally 
conjectured as the physical mechanism that induces the 
mutual electron correlations present in NSDI [9|. For 
the ellipticity value used in this paper, the field along 
the y direction drives the emitted electrons transversely 
and effectively eliminates the possibility of recollision, so 
all double ionization obtained can be regarded as orig- 
inating from sequential processes. Therefore the above 
argument for the first ionized electron also applies to the 
second ionized electron. The momentum of a resultant 
doubly charged ion equals the sum of the momenta of the 
two ionized electrons 

Px = Plx + P2x 

= -^(Ey{h) + Eyit2))^0; (9) 

Py = Ply + P2y 

= -{E,{h)+E^{t2)) = -{±Ei±E2). (10) 

The ion momentum distribution would also be ex- 
pected to have a single peak structure centered at zero 
along the x direction. What is interesting is the momen- 
tum distribution along the y direction. As Eq. ljlOp shows, 
the ion momentum distribution along the y direction is 
expected to have four peaks positioned at ±£{Ei+E2)/ijJ 
(two outer peaks) and at ±£{—Ei + E2)/oJ (two inner 
peaks) [20|,[2lj. As in the case of single ionization, we also 
get a clear relation between the ion momentum distribu- 
tion and the ionization fields of the emitted electrons. 

If the doubly charged ion momentum distribution is 
measured experimentally using COLTRIMS, projection 
onto the y direction should give four peaks, as demon- 
strated in the middle panel of Fig. [51 From the positions 
of the four peaks, noted as ±Pout and ±Pi„ with Pout > 



and Pin > 0, one gets the ionization fields of both emitted 
electrons as 



El — — [Pout ^ Pin)', 



El 



2e 



{Po. 



-* ?> 



(11) 
(12) 



The ionization times of the two electrons can also be ob- 
tained numerically using a similar method as illustrated 
in Fig. H 

The same strategy can also be applied to triple ioniza- 
tion. The ionization fields and the ionization times of all 
three ionized electrons can be obtained from the experi- 
mentally measured ion momentum distribution. Ideally, 
one would expect eight peaks along the y direction, as 
demonstrated in the right panel of Fig. [21 The positions 
of the eight peaks are labeled as ±Pi, ±i^2, i-Ps, and 
±P4, where 



(13) 
(14) 
(15) 
(16) 



There are four equations with only three unknown vari- 
ables, so knowing any three of {Pi, P2j -P3, ^4} the fourth 
one can be calculated. The ionization fields Ei, E2 and 
E'i can be deduced from the positions of the peaks 



Pi -- 


= -{-EI-E2 + E3); 
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P2 -- 


= - {+Ei - E2 + E3) ; 
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= - i+Ei + E2 + E3) . 

UJ 



Ex 
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(P2-P1); 



i^. = £(P4-P.) = £(P3 



Pi); 



P3 



4£ 



(Pi + P2 + P3 + P4). 



(17) 
(18) 
(19) 



III. NUMERICAL EXPERIMENTS 



To test the above Simpleman analytical theory, numer- 
ical experiments have been performed using the classical 
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FIG. 4: Projection of a classical ensemble onto the X1-X2 
plane. Each black dot represents a model atom. The but- 
terfly shape is a manifestation of mutual electron repulsion. 



ensemble method. This method has been described in 
detail elsewhere jl^ • 

The first step is to generate a microcanonical ensemble 
of classically modeled atoms [25]. In this paper, a model 
atom is generated with three active electrons. So far 
as we know, full-dimensional atoms with three or more 
active electrons have only been modeled using classical 
approaches [15, 26] . The ensemble is generated such that 
the total energy is fixed for each ensemble member (i.e., 
each model atom). This total energy Efot can be ex- 
pressed as 



Et. 



E ? 



VrfT^ 



E 



1 



62 



(20) 



where pi and rt are the momentum and position of the ith 
electron, and r^ is the distance between the ith and the 
jth. Note that the Coulomb potential has been softened 
with parameter a (between ion and electrons) and pa- 
rameter b (between electrons) to stabilize the classically 
modeled atom [27]. We have set a to be 1.0 a.u., to pre- 
vent autoionization, and b to be 0.1 a.u., to avoid numer- 
ical singularities. Etot is set by summing the first three 
ionization potentials, and we take -4.63 a.u., by choosing 
to model three electrons in neon. Given the total energy, 
the positions and momenta of the three electrons within 
an atom are randomly assigned. A projection of the en- 
semble phase space onto the X1-X2 plane is shown in Fig. 
m Each black dot represents a model atom (an ensemble 
member). The distribution shows a butterfly shape with 
lower probabilities in the first and third quadrants than 
in the second and fourth quadrants, a manifestation of 
mutual electron repulsion. 

Then a laser pulse is turned on and the motion of the 
electrons is governed by the TDNE 



df* 
dt 



dH 
dp' 



dp 
dt 



dH 
dr 



(21) 



The Hamiltonian including the time-dependent laser in- 
teraction is 



H = H{t) - Pt, 



3 

E 



[x^E,{t)+y,Ey{t)]. (22) 



The TDNEs are integrated numerically from the be- 
ginning to the end of the pulse. The laser field is given 
the common experimental wavelength of 780nm (oj = 
0.0584 a.u.) and an ellipticity of 0.5. The pulse has a 
sine-squared shape with full duration of 10 optical cycles 
(FWHM = 5 cycles), as shown in Fig. [3] The positions 
and momenta of electrons are recorded at each time step. 
We have defined ionization as complete when an electron 
reaches a distance of 6 a.u. from the ion core [20] and 
we have checked that a slight difference in this definition 
will not affect our discussion here. 

At each time step, the ionization criterion is applied 
to check each electron's status. If at some time step an 






FIG. 5: Left: 2D momentum distribution of singly cliarged ions at 1 PW/cm^. Center: Projection of the 2D momentum 
distribution on the x direction. Right: Projection of the 2D momentum distribution on the y direction. 



electron is detected to reach the 6 a.u. shell (note that 
an electron cannot reach this distance in the absence of 
the laser field), this time step is labeled as ii and the 
laser field strength at this time is labeled as Ei so 



Ei = JE4tiy + Ey{ti)^. 



(23) 



Due to the field ellipticity, the possibility of recoUision 
can be safely neglected (recall Fig. [T|). If at some later 
time step, a second electron is detected to reach the 6 a.u. 
shell, this time step is labeled as ^2 and the laser field 
strength at this time is labeled as E2, which is defined 
similarly as Eq. (|23p . The same strategy can also be 
appfied to the third ionized electron, and the ionization 
time and the ionization field will be labeled as t^ and 
E3 . At the end of the pulse, depending on the ionization 
results, this model atom will be classified into one of the 
four possible outcome categories: no ionization, single 
ionization, double ionization, and triple ionization. Each 
category will then be analyzed separately. 

Our numerical experiment can be "more than" a real 
lab experiment. The former can get what the latter 
can, namely, the end-of-pulse ion momentum distribu- 
tions, and it can also get what the latter cannot, namely, 
the ionization fields and the ionization times of electrons 
recorded during the pulse. Recall that the strategy of our 
Simpleman analytical theory introduced in the previous 
section is exactly to find this kind of ionization informa- 
tion from experimentally measured ion momentum dis- 
tributions. Therefore numerical experiments are ideal to 
test the validity and precision of the analytical theory: 
We start from the "experimental results" (the numeri- 
cal end-of-pulse ion momentum distributions) , apply the 
analytical theory and get the ionization fields and the 
ionization times, and then compare these values with the 
actual values (the numerically recorded values). 



IV. MULTI-IONIZATION RESULTS 

To begin. Fig. [5] shows the end-of-pulse momentum 
distributions of singly charged ions for pulses with peak 
intensity of 1 PW/cm^. One can find, consistent with 
the analytical theory, that the momentum distribution 
along the x direction peaks at zero. The appearance of 



the shoulders is due to the fact that the 2D ring shape 
shown in the left panel is not uniformly distributed. In- 
stead, it is denser in the first and the third quadrants 
than in the second and the fourth quadrants. A clear 
double peak structure can be seen along the y direction. 
The peaks are located at ±0.675 a.u. Using Eq. ([SJ of 
the Simpleman theory, the ionization field of the elec- 
tron is expected to be 0.079 a.u. On the other hand, 
as mentioned above, the classical numerical experiment 
allows one to know at what time and under what laser 
field strength an electron was actually ionized. The aver- 
aged ionization field of all single ionization events is 0.073 
a.u., which deviates less than 10% from the value inferred 
from the transverse momentum distribution. The differ- 
ence may be due to the Coulomb attraction from the 
ion core, which was taken into account all the time in 
the numerical experiment but was neglected in the an- 
alytical theory. The ionization time found numerically 
using the method illustrated in Fig. [3]is 2.76 cycles from 
the beginning of the pulse, almost exactly the same as 
the value recorded by the numerical experiment, which 
is 2.77 cycles from the beginning of the pulse. One can 
see that under our extension to an elliptically polarized 
field the Simpleman analytical theory remains valid and 
precise for electron ionization dynamics. 

The single ionization double peak structure has been 
used by Arissian, et al., to obtain the ionization field 
|28| . Circular polarization was claimed although in prac- 
tice, with such high intensity and short pulse duration, 
elliptical contaminants are usually difficult to avoid. The 
angular distribution of the 2D momentum distribution 
has been used by Eckle, et al., to measure the time that 
an electron needs to tunnel through a Coulomb barrier 

Figure [S] shows the end-of-pulse momentum distribu- 
tion of doubly charged model neon ions for peak inten- 
sity 3 PW/cm^. The 2D momentum distribution shows 
a four-band structure, corresponding to the four peaks 
when projecting onto the y axis. The momentum distri- 
bution along the x direction shows a broad single peak 
structure centered at zero. The SDI four-peak structure 
has been observed in experiment [ll|, [2l| and explained 
in detail by our classical ensemble method |2G |. 

Table |T] compares the ionization fields and the ioniza- 
tion times of the two ionized electrons, by locating the 
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FIG. 6: Left: 2D momentum distribution of doubly charged ions at 3 PW/cm^. Center: Projection of the 2D momentum 
distribution on the x direction. Right: Projection of the 2D momentum distribution on the y direction. 





Analytical Theory 


Numerical Experiment 


El 


0.093 a.u. 


0.091 a.u. 


E2 


0.25 a.u. 


0.22 a.u. 


ii 


2.27 cycles 


2.23 cycles 


i2 


5.00 cycles 


5.14 cycles 



TABLE I: Comparison of the ionization fields and the ion- 
ization times of the two electrons, from the analytical theory 
and from the numerical experiment. 



peak positions and using the analytical theory, and by 
records of the numerical experiment. One sees that the 
analytical theory fits the numerical experiment pretty 
well. The small discrepancy on ^2 between the theory 
and the numerical experiment is due to the fact that ^2 
is close to the top of the envelope, which is flat, so the 
ionization time can have a relatively large uncertainty. 

Recently, Pfeiffer, et al., have used elliptical polariza- 
tion (with ellipticity value around 0.8) to measure the 
ionization times of the two electrons in SDI [ll|. The 
method used therein is different from our method de- 
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FIG. 7: A typical triple ionization trajectory. Distances of 
the three electrons from the ion core are plotted as a function 
of time, in laser cycles. The inserted figure shows the same 
trajectory, but in a much larger space scale. 



scribed above. In [Tl| , the momenta of the two electrons 
are measured in addition to the momentum of the ion and 
the three particles are collected in coincidence, meaning 
that one must carefully check whether the two electrons 
and the doubly charged ion are actually from the same 
atom. A coincidence experiment requires a strict vacuum 
condition and a very low double ionization rate, such that 
each time only one atom is ionized, to eliminate possible 
contaminants from the ionization of neighboring atoms. 
Even when the experiment is performed with extreme 
care, false coincidences cannot be fully eliminated [301. 
In contrast, we propose in this paper that by just mea- 
suring the ion momentum, the same information of the 
ionization times of the two electrons can be obtained. No 
electron momenta are needed. No coincidence measure- 
ments are needed. Therefore the efficiency of the exper- 
iment may be substantially improved and one no longer 
needs to worry about false coincidences. 

Finally, we predict similar features in triple ionization. 
A typical triple ionization trajectory is shown in Fig. [T] 
Before emission, electrons are subject to random and fast 
collisions from the other electrons and the time scale of 
these internal collisions is much shorter than an optical 
cycle. After emission, electrons are driven away by the 
laser field and may travel for a distance on the order 
of 1,000 a.u. through the end of a typical pulse. The 
transverse momentum distribution of triply charged ions 
generated by 30 PW/cm^ pulses is shown in Fig. HI Eight 
peaks can still be recognized although the separations 
between two neighboring peaks are small. We have also 
compared the ionization fields and the ionization times of 
the three ionized electrons obtained from the analytical 
theory and the numerical experiment, as tabulated in 
table HH 



V. DISCUSSION AND SUMMARY 

In summary, we have focused on elliptical polarization, 
with which recollisions rarely happen [31], although most 
attention in strong-field atomic physics has been paid 
to recoUision-based physical processes. We show that 
elliptical polarization has the ability to reveal ionization 
information that is unreachable with linear polarization. 
Examples include the ionization fields and the ionization 
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FIG. 8: Momentum distribution of triply charged ions along 
the y direction. Eight peaks can be recognized although sep- 
arations between each two-peak pair is small. Laser peak 
intensity is 30 PW/cm^. 





Analytical Theory 


Numerical Experiment 


El 


0.070 a.u. 


0.081 a.u. 


E2 


0.34 a.u. 


0.32 a.u. 


Ez 


0.68 a.u. 


0.68 a.u. 


ti 


1.11 cycles 


1.10 cycles 


t2 


2.27 cycles 


2.23 cycles 


U 


3.76 cycles 


3.79 cycles 



TABLE II: Comparison of the ionization fields and the ioniza- 
tion times of the three ionized electrons, from the analytical 
theory and from the numerical experiment. 



times of emitted electrons. 

We have extended the Simpleman theory to include 
elliptical polarization and to predict, for the first time, 
links between the ion momentum distributions and the 



ionization fields and the ionization times of the emitted 
electrons. The ion momentum distribution along the mi- 
nor polarization direction, which is not available with 
linear polarization, is shown to contain previously un- 
expected rich ionization information and is predicted to 
have specific peak structures: a double-peak structure 
for singly charged ions, a four-peak structure for dou- 
bly charged ions, and an eight-peak structure for triply 
charged ions. The triple ionization eight-peak structure 
has not yet been observed in experiment. 

One should note that the separation between peaks 
can be enlarged by using a longer wavelength. As pre- 
dicted by Eq. ([7]) and Eq. ([T(I|. the positions of peaks 
are inversely proportional to w, thus proportional to the 
wavelength A. Longer wavelengths will be able to resolve 
close peaks that are not able to be resolved with 800nm, 
such as the ones for Ar in 21]. Substantial technical 
advancements have been made in the direction of longer 
wavelengths [32| and greater potential in exploring new 
ionization dynamics is to be expected. 

The extension of the Simpleman theory to include el- 
liptical polarization has been tested by our numerical ex- 
periments using the classical ensemble method, including 
three active electrons in the model atom. The numeri- 
cal experiments confirm the validity and accuracy of the 
Simpleman analytical theory. 
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